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Abstract Periodic cell arrangements of tracheids in Cryp- 
tomeria japonica and Agathis sp., and that of wood fibers in 
Magnolia obovata were examined by Fourier transform 
image analysis (FTIA). The angular distribution functions 
on power spectral patterns (PSPs) transformed from the dot 
maps were traced. As well as this conventional method, we 
developed the line convolution method and domain trans¬ 
formation for more correct analysis of radial and tangential 
arrangements. Cell arrangements became clear in Cryptom- 
eria and Agathis. Fluctuations were expressed by the stan¬ 
dard deviation (o) or the relative standard deviation (o r ) so 
that we could compare the characteristics of both woods. In 
Cryptomeria, the most provable distance between tracheids 
along radial files was 36.7 pm with 17% fluctuation, whereas 
the tangential interval was 28.3pm with o r of 24%. In 
Agathis, the radial diameter was 35.1 pm with o r of 23% and 
the tangential interval was 41.5 pm with o r of 23%. In Mag¬ 
nolia, the maximal periodicity of wood fibers was 15.8 pm 
and 20.5 pm along and between radial files, respectively, 
although the fluctuations could not be estimated because of 
the remarkable interruptions by numerous vessels. 

Key words Tracheid arrangement • Fiber arrangement • 
Periodicity ■ Cell arrangement fluctuation ■ Fourier 
transform 


Introduction 

Wood is an assembly of cells. When regarding xylem 
ground tissues, tracheids, and wood fibers as “structure 
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elements,” the wood structure on the transverse section can 
be understood in terms of three kinds of “shape factors,” 
namely, cell arrangements, cell shapes, and cell wall thick¬ 
ness. 1 Because the shape and arrangement are closely 
related, actual wood characteristics will become clear by 
analyzing the cell assembly pattern in addition to each cell 
property. We examine cell arrangement in this report, and 
cell shape in the following report. Cell wall thickness was 
also analyzed. 2,3 If the analysis is limited to only the struc¬ 
tural characteristics of each cell, the usual ordinary accumu¬ 
lation image analysis (OAIA) is sufficient. On the other 
hand, Fujita et al. 4 6 and Hata et al. 7 analyzed the group of 
cells as a periodic structure rather than each individual cell 
and they introduced a new technique to analyze the cell 
group using the two-dimensional Fourier transform 
method. The dot map (DM) method was developed, which 
placed a uniform point in the center of each cell on the 
transverse section to clarify the cell arrangement. 4,6 The 
most frequent (abundant) cell arrangement in the two- 
dimensional analysis was measured from the characteristic 
of the power spectrum pattern (PSP) which was output by 
fast Fourier transform (FFT)-processing of the DM. It was 
by analysis of PSP of the DM (DM-PSP) that the oblique 
arrangement of tracheids in softwood was detected for the 
first time 4 and new properties of vessel distribution were 
shown to be different from observations made with the 
traditional microscope. 8 More precise analysis of the power 
distribution demanded comparison of many specimens. 
Maekawa et al. 9 and Fujita et al. 10 developed a program 
for PSP polar coordinate analysis and demonstrated the 
tracheid and wood fiber cell arrangement by analyzing the 
power concentration, namely, the relative height of power. 
By surveying the appearance of peaks and their heights that 
reflect the cell arrangements, comparison among specimens 
became possible. However, examination of the diffusion of 
the peak was not performed. 

Before using FFT, pictures need to be transformed into 
discrete data. There are various equipment problems, in¬ 
cluding narrowness of the power brightness range. In a 
previous report, we examined the occurrence of errors due 
to picture sampling and the discrete images 1 and problems 
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Table 1 . Definitions of abbreviations and notations 

Abbreviation Definition or explanation 

or descriptor 


OAIA 

FTIA 

Structure elements 
Shape factor 
Pre-FFT 

Post-FFT 
Assembled tissue 
Dispersed tissue 

P(<t>) 

P(k) 

sNM 

^SPN 

P eq (£) and P m „(k) 
Period (x) 

Max-o 

Max-s 


0max 

k 

,v max 

tmax 

a r 

DcDM 

LcDM 


Ordinary accumulation image analysis 

Fourier transform image analysis 

Objects such as tracheids or vessels 

Shape such as diameter or area of structural elements. 

Procedures on the image domain, such as extraction of structure elements and shape factors, and also convolution to 
enhance the power level 

Procedures on the frequency domain such as the polar coordinate analysis 

Xylem tissue assembled with uniform size and shape of cells, such as tracheids in softwood 

Xylem cells such as vessels dispersed in the ground tissue 

Angular distribution function of PSP 

Frequency distribution function of PSP 

Selective net map which consist of only objective structure element, such as a body-NM and a tip-NM made from 
contours of only the cell body part or cell tip, respectively 
Shortage error of periodic number 

Frequency distribution functions of PSP along the equator and the meridian 
Periodicity distribution function 

Original maximum; the most frequent period directly calculated from the original frequency distribution function 
Peak-separated maximum; the most frequent period after the peak separation of the periodicity distribution 
function, equal to x max 
The angle that has maximum value of P (tp) 

The frequency pixel that has maximum value of P ( k ) except 0 (the origin) 

The period that has maximum value of main peak of peak-separated Period (x) 

Relative standard deviation 
Disc-convoluted dot map 
Line-convoluted dot map 


associated with Fourier transform image analysis (FTIA). 11 
It became clear that the concept of convolution on the 
images and the Radon transformation, and the 8 function, 
normally considered an error, could be used to our advan¬ 
tage. 11 Therefore, we analyzed “fluctuation,” which is not 
clearly defined, in addition to the most frequent periodicity. 
Because wood is a natural product, the arrangement of cells 
exists between periodicity and randomness. Then, for the 
present, considering the disorder from the complete period¬ 
icity as fluctuation, we evaluated the fluctuation by use of 
the power diffusion on PSPs. 

The terms and abbreviations associated with OAIA and 
FTIA were announced in previous reports 1-3,11 and parts of 
them are shown here. Table 1 shows a list of abbreviations 
used throughout this report. 


Materials and methods 

Specimens and microscopic preparations 

For verification of the results, we selected uniform samples, 
namely, the early wood of Cryptomeria representing little 
change in cell arrangement and shape, and Agathis, a tropi¬ 
cal softwood with a relatively constant structure and no 
annual rings. On the other hand, Magnolia was chosen as an 
example of hardwood fiber, because its cell arrangement 
was the most disperse among many Japanese hardwoods 
tested by Maekawa et al. 9 and Fujita et al. 10 The microscopic 
preparations were described previously. 111 The process in¬ 
volved treatment of the wood block with 55% FIF, epoxy 


resin embedding, slicing of 5-pm-thick transverse sections, 
safranine staining, and balsam mounting. 

Original picture sampling and equipment 

When quantitatively analyzing the range of fluctuation, it is 
necessary to obtain a large number of structure elements. 
We used very low-magnification photographic conditions of 
2X (objective lens) and 2.5X (projection lens) for Cryptom¬ 
eria and Agathis, and 4X and 2.5X for Magnolia , to sample 
as many cells as possible at a time. The photographic 
pictures were transformed to digital images by a high- 
resolution film scanner at 1280 dpi. 1-12 The limitation factor 
of the resolutions was R a (4.0 pm in Cryptomeria and 
Agathis, and 2.0 pm in Magnolia). We selected Luzex III 
(Nireco) as the main image processor, because it expresses 
the Perceval equation without suppression of power bright¬ 
ness. 11 The input method of the transverse section image 
was the same as previously reported. For example, the di¬ 
rection of radial growth was adjusted to coincide with the 
x-axis on the display screen. 1-11,12 

Extractions of structure element and shape factor, and 
some operations 

We made input pictures into discrete images, namely, cell 
wall maps (WM) or cell lumen maps (LM). In Cryptomeria 
and Agathis, a small number of ray parenchyma cells 
was involved, although almost all cell lumens were those 
of tracheids. Because ray lumens were so slender along the 
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Fig. la,b. Result of Fourier transform image analysis (FTIA) proce¬ 
dure on a cross sectional view of Cryptomerici japonica. a A part of the 
discrete image of the original picture input to the image analyzer, b 
Discrete image of dot map-power spectrum pattern (DM-PSP), show- 


x-axis, they were excluded by size filters. Axial parenchyma 
cells were ignored. The example of Cryptomeria is shown in 
this article (Fig. 1), while those of Agathis, shown in previ¬ 
ous reports 111 are omitted here. In Magnolia, lumens were 
extracted and those of vessels and rays were deleted to give 
only wood fibers. 

For extraction of the shape factor, each lumen was sim¬ 
plified to the central point of cells that constituted one pixel, 
namely, the dot map (DM) according to the method pro¬ 
posed by Maekawa et al., 9 and it was cut off in a circular 
window. 11 To reinforce power, dots were expanded uni¬ 
formly using the method of Maekawa et al. 9 This map and 
the PSP are called the disc-convoluted dot map: DM, DM- 
PSP, respectively. We used the polar coordinate analysis 
program developed by Maekawa et al. 9 and Fujita et al. 10 for 
analysis of PSP. In addition, we developed a new convolu¬ 
tion and analysis system as follows. 


Results and discussion 

Power reinforcement and concentration effects by line 
convolution 

Dot maps with a large fluctuation in their arrangement 
resulted in only a weak power level, even if there were more 
than 1000 dots. Until now, the DM method has been used for 
power reinforcement with pre-FFT (see Table 1). An ex¬ 
ample for Cryptomeria is shown in Fig. 1. This method is 
convenient for visual judgment of the PSP because the PSP 
is output onto an Airy disc preserving the two-dimensional 
information. However, the power decreases with the fre¬ 
quency pixels ( k ) and it generates the so-called missing 
order of an Airy disc around 250, 11 so quantitative analysis of 
a frequency distribution function will be limited to the low 
frequency domain. Although Maekawa et al. 9 performed a 
“deconvolution operation” using a random dot map PSP for 
the offset of this weak point, this method was very complex. 
Therefore, we developed a very simple and effective tech¬ 


ing the power along the meridian on which spots can be recognized 
(arrows). Four spots ( arrowheads) are located on the two striations and 
suggest a honeycomb assembly of tracheids. The lower half was output 
at lower brightness than for the upper half 


nique based on the principle of the Radon transformation. 
In the example shown in Fig. 2, each dot was extended up 
and down for four pixels (nine pixel lengths in total) and the 
power was amplified nine times by convolution. The opera¬ 
tion by which dots are extended vertically or horizontally is 
very easy for the image processor. Also, another effect can 
be added on the power distribution that is spread into the 
two dimensions in case of DcDM-PSP. That is PSP of line 
convoluted dot maps (LcDMs) is concentrated and 
amplified near the equator line by 1024/9 (Fig. 2b). As these 
extended lines are one pixel in width, periodic information is 
maintained in the x-direction by the Radon transform effect. 
That is, each extended line has a strong point as the 6 
function as described in the previous report. 11 This convolu¬ 
tion method was so useful that LcDM-PSPs were used to 
examine the fluctuation of arrangement in radial and tan¬ 
gential directions. 

On the other hand, we contrived another technique to 
extract only the wood fiber bodies and enhance the power 
(Fig. 3). In the usual net map (NM), the cell boundary is 
extracted. In softwoods, each neighboring tracheid shares 
this boundary line. However, when adjacent cells are of 
different types, the cell boundary line must be made clear. 
There are many vessels and ray parenchyma cells in hard¬ 
woods. Moreover, there is a difference between the body 
part and the tip part. Therefore, we developed the tech¬ 
nique described below to determine the origin of body 
parts. At first, the NM was expanded to five pixels wide (this 
width can be optionally selected according to software of 
the equipment) and the reverse image was made. Size filters 
removed vessels and rays, and then the body part (Fig. 3b) 
of large size was divided from the small tip part with triangle 
or diamond form by a size filter. The ratio of the tip part to 
the body part was about 0.3. When these lumens are pro¬ 
cessed to dots, the selective dot maps are called “body-DM” 
and “tip-DM,” respectively. On the contrary, the differen¬ 
tial filter can extract the edge of lumens, although the posi¬ 
tions are shifted from the true wall-lumen boundaries. The 
main purpose of this pre-FFT processing is cell shape analy¬ 
sis, but the reinforcement effect of the power was remark- 
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Fig. 2a-d. Line convolution ef¬ 
fect in Agathis sp. a Line convo¬ 
luted dot map (LcDM) consisted 
of dots elongated vertically to 
nine pixels. Lines are dilated by 
two pixels in the horizontal direc¬ 
tion for clarity, b LcDM-PSP. c 
Frequency distribution function 
on the equator of LcDM-PSP 
[l\ q (k) = P- ss^ 4 {k)\. d Fre¬ 
quency distribution function on 
the meridian of LcDM-PSP 
[Pmer(F) = ] ‘ Note 

LcDM-PSP in b shows that the 
power is concentrated around the 
equator. Although large peaks 
are clear in c and d, additional 
peaks are also detected slightly 
( arrows ) 



O Frequency distribution function of Agathis sp.’s 
LcDM-PSP P. 5S4S< (A) 



b 



Sf: 32pix. 

mm - 


Frequency distribution function of Agathis sp.’s 
LcDM-PSP P, SS4£m (A) 



Frequency pixels [pix.] 


able. This is called a “body-NM,” which is one of the selec¬ 
tive NM (sNM). The power is expressed around k — 200 
corresponding to paired boundaries. 

Visual analysis of DcDM-PSP 

Interpretation of PSPs that are presented on the two- 
dimensional polar coordinates of the frequency domain 
needs practice because phase information on their original 
image domain was erased. Moreover, size information is 
expressed by a reciprocal. 11 However, when an input picture 
is simplified to DcDM by pre-FFT, the visual understanding 
of PSP is not difficult (Fig. 1). Tracheids and wood fibers 
produced from the fusiform initials by the periclinal division 
make a series of cells, that is, radial file. The arrangements 
of intraradial files appear on the equator line of DM-PSP, 
whereas those of interradial files are expressed on the 
meridian. 11 

In Cryptomeria , two types of striations were observed; 
one was the meridian striation, the other was a paired stria- 
tion at frequency pixel k = 118 stretched in parallel to the 
former (Fig. lb). There were two concentrated parts of 
power at k = 150 on the meridian striation (arrows in Fig. 
lb). The power spots were also detected at two points on 
paired ones (arrowheads in), respectively, then six center¬ 
ing parts in the amount existed. This result was the same in 
Agathis . 11 It suggested that tracheids were put together like 
a honeycomb ambiguously, but the arrangement surely 
existed. 4-6 

In case of Magnolia , the power having a ring was ob¬ 
served as in an analysis example in a previous report by 
Maekawa et al. 9 However, when observing in detail, stria¬ 
tions that stretched along the direction of the meridian were 
observed. Figure 3d shows DcDM-PSP that was generated 


from the body-DM (Fig. 3c) and in this PSP the striation is 
clearly observed (arrows). On the other hand, the same 
characteristics were also observed in the tip part PSP (PSP 
was omitted). These facts suggest that the body part and tip 
part together were arranged in the radial direction. No 
other peculiar power was discovered from the tip part PSP. 
The power concentrations were also observed in the ob¬ 
lique direction, in which arrowheads pointed in the body 
part PSP shown in Fig. 3d. However, the power was so 
diffuse that specifying a maximal part was not easy. Because 
the power was weak in DcDM-PSP, we examined it more 
in the sNM-PSP (Fig. 3f) generated by the FFT process 
from the sNM, namely, the body-NM shown in Fig. 3e. A 
power distribution similar to the case of DcDM-PSP was 
observed in the low-frequency domain (arrowheads) and 
moreover the power like a ring appeared at about k = 200 
in the middle-frequency domain. Other spots in the equator 
direction were detected in addition to the spots on the 
meridian and the oblique direction. Eight spots could be 
specified (arrows) due to the effect that boundary lines 
composed of two lines with a five-pixel interval reinforced 
and concentrated the power in a ring as described. 

Reconstruction of the most frequent model of tracheid 
and wood fiber arrangements 

Tracheids and wood fibers produced from the fusiform ini¬ 
tials by the periclinal division make a radial file. The cell 
arrangement of the intraradial file appears on the equator 
of the DM-PSP. For more precise quantitative analysis and 
the concrete reconstruction of cell arrangements, we ob¬ 
tained the frequency distribution function P ( k , <j>) of 
LcDM-PSPs from the sector —5 SI i) § 4 including the 
equator, and 85 S 0 S 94 including the meridian (Fig. 2c,d). 
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Fig. 3a-f. The general procedure of FTIA to extract only the body part 
of wood fibers in Magnolia obovata. a Part of discrete image of original 
picture, b Lumen map of body part extracted after the removal of fiber 
tips, vessels, and rays by size filters, c Disc convoluted dot map 
(DcDM) from b. d DcDM-PSP. e Selective net map (sNM), especially 
called a body-NM, produced by differentiation of lumens in b. f sNM- 
PSP. DcDM-PSP shows that the power is distributed circularly and 
some spots are discerned on the circle (arrows and arrowheads). A 
lower half was outputted at lower brightness than for the upper half in 
d. sNM-PSP shows similar spots (arrowheads) to d and further spots on 
a higher frequency power ring ( arrows) 


As described in a previous report, 11 we designated the 
brightness distribution function of the map as B (x, y ), the 
Fourier spectrum as Four.(k, </>), and the power spectrum as 
P ( k , <j>). Flere, we propose the periodicity distribution func¬ 
tion Period (x, <j>) of the map. Because discussion is often 
focused on a certain direction, we express Four.(k, 0), P (k, 
(j>), and Period (t, <j>) as Pour.(k). P ( k ), and Period (x), 
respectively. Analysis starts from P ( k) (Fig. 2c). P (k) was 
changed to Four.(k) by calculating the square root of the 
vertical axis (Fig. 4a), and changed to Period (x), which was 
obtained by translating the horizontal axis reciprocally 
(Fig. 4b), namely, x[pm] = 1024/A: [pixels] X 4.0 [pm/pixel]. 
Period (x) was submitted to the peak separation soft¬ 
ware Sigma Plot to obtain the maximum x (x max ). Separa¬ 
tion assumptions and factors are as follows: a function 
curve is separated as a normal distribution function and a 
curve consists of a main (primary) peak and a sub (high- 


a 


Frequency distribution function of Agathis sp.’s 
LcDM-PSP Four .-(*) 



0 50 100 150 200 250 300 350 400 450 500 

Frequency pixels [pix.] 


b 


Periodicity distribution function of Agathis sp.’s 
LcDM-PSP Period (r) 



0 10 20 30 40 50 60 70 80 90 100 

Period [/f m] 



0 10 20 30 40 50 60 70 80 90 100 

Periodic m] 


Fig. 4a-c. Periodicity analysis in r-direction in Agathis sp. a Fourier 
spectrum calculated from the frequency distribution function on the 
equator of LcDM-PSP (Fig. 2c). b Periodicity distribution function on 
the horizontal axis transformed from a. c Peak separation of b by using 
the Sigma Plot software package. Note that although the periodicity 
distribution function was well separated, further peaks remain (arrow) 


order) peak. Peak height (a); « nla]n > a sub , peak position (x); 
x mam > X su b , constant value (c); c maln = c sub . The correlation 
coefficient R was examined, R = 0.945 (Fig. 4c). Finally, 
from the x max cell position was calculated in the radial and 
tangential directions. In Table 2, Max-s (peak-separation 
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Table 2. Values on tracheid and fiber arrangements calculated from the angular and frequency distribution functions of DM-PSP 


Samples 

Power spectral pattern 


Calculated periodicities 


^"(degrees) 

k ma x c (pixels) 

Max-s d ± o(pm) 

a/(%) 

Max-o f ( pm) 

Cryptomeria japonica 

0 (r-direction) 

118 

36.7 ± 6.1 

17 

34.7 

if = 6000 

90 (t-direction) 

150 

28.3 ± 6.9 

24 

27.3 


156 (oblique direction) 

127 

- 

- 

32.3 

Agathis sp. 

0 (r-direction) 

123 

35.1 ± 8.0 

23 

33.3 

n = 6600 

90 (t-direction) 

98 

41.5 ± 9.5 

23 

41.8 


159 (oblique direction) 

130 

- 

- 

31.5 

Magnolia obovata 

0 (r-direction) 

130 

- 

- 

15.8 

n = 3100 

90 (t-direction) 

100 

- 

- 

20.5 


46 (oblique direction) 

86 

- 

- 

23.8 


“Number of structure elements 

b The angle that has maximum value of angular distribution function 

'The frequency pixel that has maximum value of frequency distribution function except 0 (the origin) 

d The most frequent period after the peak separation of the periodicity distribution function, equal to T max . The curve was separated as a normal 
distribution and o was calculated mathematically 
'Relative standard deviation (100 o/Max-s) 

f The most frequent period directly calculated from the original frequency distribution function 


maximum) denotes x max . A concentration of power in the 
oblique direction was detected in both softwoods and hard¬ 
woods. 9,10 ' 13 This indicated that the cell was arranged in the 
oblique direction, and we measured the maximum angle 
(0max) from the angular distribution function P (0) of 
DcDM-PSPs (Table 2). 

In the same procedures conducted with Magnolia, k max , 
which is the frequency pixel that has maximum value of 
P ( k ), was obtained from the body-LcDM-PSP (Table 2). 
However, the peaks on both the equator and the meridian 
were different those for softwood. They were prominent on 
only one part of them (graphs not shown). k max in Table 2 
shows this prominent part and Max-o (original maximum) is 
the period that is directly calculated from k max of the PSP by 
a conventional technique. The arrangement of the oblique 
direction was obtained from 0 max , which was weak but was 
detected in P((j)) of DcDM-PSP. Although this angle could 
not be detected in the visual judgment of DcDM-PSP 
(Fig. 3d), concentrated powers in the oblique direction 
could be distinguished in the visual judgment of sNM-PSP 
(arrowheads in Fig. 3f). The angles of these powers were 
almost the same as 0 max of DcDM-PSP. 

A model of cell arrangement with the highest possibility 
was reconstructed from the information in Table 2 and is 
shown in Fig. 5. A model restored by the central point of 
cells is called a dot model. However, the dot model is a very 
different image from that of the actual transverse section, 
which prompted us to prepare the block model. In soft¬ 
wood, tracheids were considered an assembled tissue, be¬ 
cause cell walls were situated on the middle of the dots, and 
a line was used to divide each dot at the center (Fig. 5a,b). 
Cell wall angles cannot be determined until angle analysis is 
done (this will appear in the following report). Tracheid 
shape was supposed to be rectangular considering the radial 
files. This model is called a block model because it re¬ 
sembles a brick wall. Seeing block models of Cryptomeria 
and Agathis, they are rectangular. When the comparison is 
made between sideways and longways, the former is longer 
sideways and the latter is shorter sideways, respectively. 


Thus, characteristics among species can be clarified. On the 
other hand, the dot model of Magnolia was reconstructed 
only with the body part of wood fiber, because the tip part 
was removed. Therefore, a block model was not recon¬ 
structed in this step. 

Experimental conditions for fluctuation analysis 

Dot models and block models were restored from the maxi¬ 
mal part of the power on DcDM-PSP, LcDM-PSP, and so 
on. These models are thought to represent the most fre¬ 
quent cell arrangements extracted from the whole sampling 
area. However, actual PSPs show high diffusion, being dif¬ 
ferent from the case of X-ray diffraction of crystals. This 
means that actual cell arrangements in woods are not uni¬ 
form. They have considerable “fluctuation” as is usual with 
natural products. We are sure the fluctuation is related to 
the species and age of trees, and also affects wood proper¬ 
ties. Maekawa et al. 9 obtained the angular distribution func¬ 
tion and the frequency distribution functions in 17 species 
and then obtained the relative peak heights to whole power 
using the deconvolution and normalization method. This 
was followed by work by Fujita et al., 10 and Diao et al. 13 
covering many specimens. This method is effective for com¬ 
parison between specimens but does not really relate to the 
essence of fluctuation. Also, most of peaks were in the low- 
frequency bandwidth; the concern over frequency errors 
/% 1 is large. Considering these points, we set the environ¬ 
ment to control the occurrence of various errors to within 
10% as previously reported. 1,11 Analysis conditions were 
chosen to minimize the frequency error (£ F ), considering 
the optical resolutions ( R 0 ) and discrete resolution 
(Ry ,). 111,12 Enormously abundant cells can make fine scan¬ 
ning possible to determine frequency distributions. Each 
cell was sampled at the size of about 10 pixels at low magni¬ 
fication and the negative influence of R v was avoided. 11 

Next, we examined diffusion of peak caused by “short¬ 
age of periodic number” (spn). When there are few 
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3 Cryptomeria japonica 



C Magnolia obovata 


• • 


20.5 1U m 
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M 

15.8/i m 



• • 

0 = 46° 


100/i m 


Fig. 5a-c. Reconstruction of cell arrangements and fluctuations, a Dot 
model and block model of Cryptomeria japonica. b Models of Agathis 
sp. c Dot model of Magnolia obovata. The most abundant period, the 
standard deviation, the relative standard deviation, and oblique ar¬ 
rangement (0) were added to models. Values are selected from max-s 
in Cryptomeria japonica and Agathis sp. and max-o in Magnolia 
obovata (see Table 2) 


numbers of periods even if the period is not disordered, 
power will inevitably diffuse due to the uncertainty prin¬ 
ciple. 14 The Scherrer formula of the X-ray diffraction is 
based on this effect, and, conversely, crystal width can be 
calculated from the diffusion of peaks. Periodic number 
depends on specimen characteristics, for example, in hard¬ 
wood, vessels interrupt fiber rows. We designated the error 
of spn (£ spn ) as a ratio of brightness value in the pixel next 
to the frequency pixel that had the maximal brightness 
value. /i SPN was ignored when the period number was 50 or 
more and fs SPN became less than 10% by use of a vertical 
line pattern as in a previous report. 11 Because the cell num¬ 
ber of this examination was about 6000 in Cryptomeria and 
Agathis and about 3000 in Magnolia (Table 2), i.e., about 80 
and 55, respectively, in one dimension, it was thought that 
power diffusion of the PSP reflected the fluctuation be¬ 
tween each dot’s distance, not the spn error. 


Detection of fluctuations in three species 
and conclusions 

Power diffusion on LcDM-PSP of Agathis , for instance, at 
P 5 S ^ S4 (k), namely, / J eq (k) is shown in Fig. 2c. The diffusion 
of power is expected to originate from the probability distri¬ 
bution of periods. There was a large peak in k max = 123 and 
more peaks were detected in about 250 and 380, which were 
multiples of 123. These were considered to be subpeaks, i.e., 
the second and third peaks in the diffraction. It must be 
noted that P (k) existed in the frequency domain, while the 
probability distribution of periods is in the image domain. 
Then, the domain transformation was introduced. P{k) was 
transformed to Period (x) as already described. Subpeaks 
were overlapped at high frequency in the case of P (k), and 
at low frequency in the case of Period (x). If excluding this 
influence of subpeaks, the form of the main peak will be 
equivalent to the frequency distribution of distance among 
dots. That is, the form of the peak represents the degree of 
fluctuation. The peaks become high and sharp when dots 
are periodically arranged, whereas the disorder of the 
period, i.e., high fluctuation, results in a wide peak. To 
evaluate peak diffusion, we used the normal distribution 
analysis that was established mathematically, and the stan¬ 
dard deviation o was selected. As a result of peak separa¬ 
tion, the main peak was separated from a secondary peak 
(Fig. 4c). The correlation coefficient R was poor because the 
higher-order (second and third) peaks were considered as 
one subpeak. Relative standard deviations o r (100 o/Max-s) 
were also calculated for comparisons among materials 
(Table 2 and Fig. 5a,b). Figure 5a,b was reconstructed from 
Max-s and Fig. 5c was constructed from Max-o. We pro¬ 
posed o r as the index of an arrangement’s fluctuation. When 
the fluctuation is evaluated and expressed as o r in this way, 
cell arrangement disorder can be compared among tree 
species or among organizations. 

Regarding R eq (k) (see arrows in Fig. 2c and Fig. 4), there 
is a small peak on the main peak. This implies two periods 
along the radial direction. When the input image was ob¬ 
served once more, some latewood-like bands of about 
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30 [im in the radial direction and with a slightly thicker wall 
were present in Agathis. After such tracheids were cut off 
before dot mapping, the small peaks disappeared. In Cryp- 
tomeria sampled only from earlywood, such a small peak 
was not detected. There was a peak in addition to the main 
peak on P mN {k) (arrow in Fig. 2d). This peak of the fre¬ 
quency domain implies a shorter period between radial 
files. In Cryptomeria, a shorter period also existed. These 
periods were perhaps caused by a decrease in the tangential 
diameter of the tracheid tip part. Max-s will be examined by 
OAIA and the results will be described in the following 
report. 

As for the oblique arrangement of cells, we mentioned 
only the most frequent direction. However, it was expected 
that the fluctuations in any direction could be expressed by 
the synthesis of both vectors because o r of radial and 
tangential directions were already detected. Future studies 
are needed to determine the most suitable direction for 
analysis. 

For Magnolia , the period of dots was vague, and power 
was wide and diffusive as mentioned in the previous re¬ 
ports. 910 When we made wood fiber images of only the body 
and the tip was excluded, spots were clearer in the DcDM- 
PSP and sNM-PSP (Fig. 3d,f). However, because it was 
dimmer than in Agathis and Cryptomeria , the maximal part 
was detected but the fluctuation could not be examined. In 
hardwoods, wood fiber arrangement was inevitably dis¬ 
turbed because of the presence of vessels. E s PN must be 
reconsidered, because there were about five periodic num¬ 
bers in Magnolia , and the main peak was diffuse. Hard¬ 
woods have many properties that are different from those 
of softwoods and analysis is more difficult. However, be¬ 
cause Magnolia has the most ambiguous arrangement of 
characteristics, we want to examine more species in greater 
detail. 

Many cells were analyzed at a time in this report and 
characteristics of cell arrangements were expressed by use 
of the most frequent periodicity and “fluctuation.” A series 
of techniques to reconstruct cell shapes will be completed 
after the direction information of the walls is examined as 
well as these factors in the following report. Various 
examinations, namely, feature extraction, feature compari¬ 
son in many species, a grasp of the shape change within 
annual rings, and so on, can be performed using this 
method. The meaning of cell shape will finally be clear. 
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